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ABSTRACT 



in 

The recently discovered subdwarf B (sdB) pulsator KIC 1 1558725 is one of the 16 pulsating sdB stars detected in the Kepler field. It 
features a rich g-mode frequency spectrum, with a few low-amplitude p-modes at short periods. This makes it a promising target for 
5^ ' a seismic study aiming to constrain the internal structure of this star, and of sdB stars in general. 

We have obtained ground-based spectroscopic radial-velocity measurements of KIC 1 1558725 based on low-resolution spectra in the 
Balmer-line region, spanning the 2010 and 2011 observing seasons. From these data we have discovered that KIC 11558725 is a 
binary with period P= 10.05 d, and that the radial-velocity amplitude of the sdB star is 58 km s _I . Consequently the companion of the 
sdB star has a minimum mass of 0.63 Mq, and is therefore most likely an unseen white dwarf. 

We analyse the near-continuous 2010-201 1 Kepler light curve to reveal the orbital Doppler-beaming effect, giving rise to light varia- 
tions at the 238 ppm level, which is consistent with the observed spectroscopic orbital radial-velocity amplitude of the subdwarf. 
(N ' We use the strongest 70 pulsation frequencies in the Kepler light curve of the subdwarf as clocks to derive a third consistent mea- 

surement of the orbital radial- velocity amplitude, from the orbital light-travel delay. The orbital radius a s( jB sin/ = 11.5^q gives rise 
to a light-travel time delay of 53.6 s, which causes aliasing and lowers the amplitudes of the shortest pulsation frequencies, unless the 
C**) ■ effect is corrected for. 

We use our high signal-to-noise average spectra to study the atmospheric parameters of the sdB star, deriving T cff = 27 910 K and 
logg = 5.41 dex, and find that carbon, nitrogen and oxygen are underabundant relative to the solar mixture. 

Furthermore, we analyse the Kepler light curve for its pulsational content and extract more than 160 significant frequencies. We 
investigate the pulsation frequencies for expected period spacings and rotational splittings. We find period-spacing sequences of 
spherical-harmonic degrees C=\ and £=2, and we associate a large fraction of the g-modes in KIC 1 1558725 with these sequences. 
From frequency splittings we conclude that the subdwarf is rotating subsynchronously with respect to the orbit. 



Key words, stars: subdwarfs - stars: early-type, binaries, oscillations - stars: variable: subdwarf-B stars - stars: individual: 
KIC1 1558725 

1 . Introduction majority of the sdB population (see Heber 2009l for a detailed 

review). 

The hot subdwarf B (sdB) stars populate an extension of the hor- 
izontal branch where the hydrogen envelope is of too low a mass Several extensive radial-velocity surveys have targeted the 
to sustain hydrogen burning. These core helium burning stars sdB stars, with thejmost recent large sample explored by 
must have suffered extensive mass loss close to the tip of the red i Copperwheat et al. | dZOlJJ). They find that -50% of all sdB stars 
giant branch in order to reach this core/envelope configuration. reslde ln short-period binary systems with the majority of corn- 
Binary interactions, either through stable Roche lobe overflow P amons bem 8 whlte dwarf ( WD ) stars - A recent compilation 
or common envelope ejection, are likely to be responsible for the , of such short^enod systems can be found in Appendix A of 

Geieretal. (2011), w hich lists a total of 89 sy stems. Adding 

* Based on observations obtained by the Kepler spacecraft, the Kitt 18 new systems from |Copperwheat et al.| (|201l|) brings the to- 

Peak Mayall Telescope, the Nordic Optical Telescope and the William tal well above a hundred sdBs with periods ranging from 0.07 

Herschel Telescope to 27.8 d. These systems are all characterised by being single- 
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Fig. 1. Mean spectra from KP4m, NOT, WHT (bottom to top), 
offset in flux for clarity. The top panel is a zoom-in of the WHT 
spectrum, demonstrating some of the stronger lines of heavy el- 
ements, such as Mgll 4481 A and Silll 4552,4567 A. 



lined binaries, i.e. only the sdB stars contribute to the optical 
flux, which directly constrains the companion to be either an 
M-dwarf or a compact stellar-mass object. Binaries with com- 
panions of type earlier than M are double-lined and also read- 
ily ident ifiable from a combinat ion of optical and infrared pho- 
tometry. |^ei&Stie^i (|2003) find that -50% of all sdB stars 
have IR excess and must have a companion no later than M2. 
Radial velocity studies targeting these double-lined stars have 
had a hard time detecting orbital periods, indicating that they 
m ust be exceedingly long. A rece nt breakthrough was made 
by I0stensen & Van WinckeH d2012l) using high-resolution spec- 
troscopy of a sample of eight bright subdwarf + main-sequence 
(MS) binaries detecting orbital periods spanning a range from 
~500 to 1200 d with velocity amplitudes between 2 and 8 km s _1 . 
The period distribution of these different types of binary systems 
are important in that they can be used to constrain a number of 
vaguely defined parameters used in binary population synthesis 
models, including the common envelope ejection efficiency, the 
mass and angular momentum loss during stable mass transfer, 
the minimum core mass f or helium i gnition, etc. The seminal 
binary population study of lHan et alJ d2002l 120031) successfully 
predicts many aspects of the sdB star population, but the key 
parameters have a wide r ange of possib l e valu es. A recent popu- 
lation synthesis study bv lClausen et alJ (120121) explores the pos- 
sible populations of sdB+MS stars and demonstrates how the 
entire population can change with different parameter sets, but 
does not deal with sdB+WD binaries. 

A theoretical prediction of the existence of pulsations in sdB 
stars, due to an opacity bump associat ed with iron ionisatio n 
in subphotospheric layers, was made by Charpin et et al.l(ll997l) . 
Since both p and g-mqde pulsations wer e discovered in sdB stars 
(Kilk enny et alJll997l : lGreen et al.ll2003l) . there has been a focus 
on the possibilities to derive the internal structure and to put con- 
straints on the lesser known stages of the evolution by means of 
asteroseismology. Currently the immediate aims of asteroseis- 
mology of sdB stars are to derive the mass of the Helium-burning 
core and the mass of the t hin Hydrogen envelope around the core 
(e.g. iRandall et alj|2006|). the rotation al frequency and internal 
rotation profile (Charpinet et al. 2008), the radius, and the com- 



position of the core (e.g. lVan Grootel et al.l2010t fCharpinet et al. 
l201lh . 

Recent observational success has been achieved from splen- 
did light curves obtained by the CoRoT and Kepler spacecrafts, 
delivering largely uninterrupted time series with unprecedented 
accuracy for sdB stars. Overviews of the Kepler survey s t age re - 
sults for sdB stars were given bv IQstensen et al.l (1201 Obi 1201 lb . 
and case studies re vealing dense pul satio nal frequency spec - 
tra are presented by iReed et al.l (1201 Ol) and iBaran et all (1201 ll) . 
From Kepler data it has become clear that t he ff-modes in sdB 
stars can be identified from period spacings (IReed et al.ll201 lb . 
Earth-size planets stripped from their outer layers have been 
found around the pulsating sdB star K IC 05807616 (KOI-55, 
KPD 1943+4058. l(Iharpinet etal1l201 ll) . the star being also the 
subject of the first seismi c study of an sdB star in the Kepler field 
(IVan Grootel etalj|2010b . 

Kepler sdB+dM binaries with pu lsating subdw ar f com - 



ponents have been presented by Kawaler et ail (|2010), 



Oste nsen et all d2010al) . and by iPabloetalJ d2011l) . White 



dwarf companions in close Kepl er binaries are presen ted by e.g. 
iBloemen etail d201 UI20T2I) and lSilvotfi et all d2012l) . 

Our target, KIC 11558725 or J19265+4930, is one of the 
16 pulsating sdB stars detected in the Kepler field. The Kepler 
magnitude of KIC 11558725 is 14.95, and the B-band magni- 
tude is about 14.6, making it the third brightest in the sample. 
A first description of the spectroscopic properties, and the pul- 
sational frequency spectrum a s found from the 26 d ay Kepler 
survey dataset, was given by IQstensen et al.l (1201 ll) . with the 
source showing frequencies in the range of 78-391 yuHz. Based 
on this relatively short data set already 36 frequencies were iden- 
tified, showin g the poten t ial of this star for a seismic study. 
Subsequently, Bara n et al.l (1201 ll) derived 53 frequencies in total 
from the Kepler survey data, and the freq uencies were iden tified 
in terms of spherical-harmonic degrees bv lReed et al.l (1201 ll) . As 
a consequence, the star was observed by Kepler from Q6 on- 
wards. At the time of writing, we have analysed data from the 
five quarters Q6 to Q10. We present the full frequency spectrum 
resulting from these 15 months of short-cadence Kepler obser- 
vations. 

In this paper we present our discovery of the binary nature of 
KIC 11558725 based on low-resolution spectroscopy. This ob- 
ject was sampled as part of a spectroscopic observing campaign 
to study the binary nature of the sdB population in the Kepler 
field, for which some prelimin ary results have already been pre- 
sented by Telting e t al.l (120121) . Part of the data from this cam- 
paign were presente d in case studies on the close sdB + WD bi- 
nary KIC 06614501 dSilvottietalJl2012h . and the p-mode pul- 
sating sdB KIC 10139564 (Baran et al. 2012). 

From our new spectra of KIC 1 1558725 we solve the orbital 
radial-velocity amplitude, and derive a lower limit of the com- 
panion of the sdB star, which is most likely an unseen white 
dwarf (Sect. 2). We use the average spectrum to study the at- 
mospheric parameters in detail. We show that the orbital Kepler 
light curve reveals strong evidence for Doppler beaming that 
results in light variations at the 238 ppm level, consistent with 
theoretical predictions, again allowing us to make an indepen- 
dent measurement of the orbital radial-velocity amplitude (Sect. 
3). We extract 166 pulsational frequencies from the Kepler light 
curve (Sect. 4), and show that the orbit has an appreciable ef- 
fect through the light-travel time on the observed phases and 
frequencies of these pulsations, which in fact allows us to make 
a third independent measurement of the orbital-radial velocity 
amplitude (Sect. 5). In the final sections of this paper we discuss 
pulsational period spacings and frequency splittings, aiming to 
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Fig. 2. Radial-velocity curve from observations with the KP4m Mayall, the Nordic Optical, and the William Herschel telescopes. 



identify the spherical-harmonic degree of the modes and to dis- 
close the rotation period of the subdwarf in KIC 1 1558725. 



2. Spectroscopic observations 

Over the 2010 and 2011 observing seasons of the Kepler field 
we obtained altogether 35 spectra of KIC 1 1558725. 

Low-resolution spectra (R w 2000-2500) have been collected 
using the Kitt Peak 4-m Mayall telescope with RC-Spec/F3KB, 
the kpc-22b grating and a 1 .5-2.0 arcsec slit, the 2.56-m Nordic 
Optical Telescope with ALFOSC, grism #16 and a 0.5 arcsec slit, 
and the 4.2-m William Herschel Telescope with ISIS, the R600B 
grating and 0.8-1.0 arcsec slit. Exposure times were 600 s at 
KP4m and WHT, and either 600 s or 300 s at the NOT The re- 
sulting resolutions based on the width of arc lines is 1 .7 A for the 
KP4m and WHT setups, and 2.2 A for the setup at the NOT. See 
TableQ]for an observing log. 

The data were homogeneously reduced and analysed. 
Standard reduction steps within IRAF include bias subtraction, 
removal of pixel-to-pixel sensitivity variations, optimal spectral 
extraction, and wavelength calibration based on arc-lamp spec- 
tra. The target spectra and the mid-exposure times were shifted 
to the barycentric frame of the solar system. The spectra were 
normalised to place the continuum at unity by comparing with a 
model spectrum for a star with similar physical parameters as we 
find for the target (see Sect. 3.2). The mean spectra from each of 
the three telescopes are presented in Fig. 1. 

2.1. Radial velocities and orbit solution 

Radial velocities were derived with the FXCOR package in 
IRAF. We used the Hy, H6, H£ and Hrj lines to determine the 
radial velocities (RVs), and used the spectral model fit (see next 



section) as a template. See TableQ]for the results, with errors in 
the radial velocities as reported by FXCOR. The errors reported 
by FXCOR are correct relative to each other, but may need scal- 
ing depending on, amongst other things, the parameter settings 
and the validity of the template as a model of the star. As our 
fit results in a ;^ 2 -value close to unity (see Table 2) we trust that 
scaling of the FXCOR errors is not necessary and that the de- 
rived errors of the fit parameters are adequate. 

Assuming a circular orbit we find an orbital pe- 
riod of 10.0545(48)d, with a radial-velocity amplitude of 
58.1(1.7)kms _1 for the subdwarf. See Table [2] for the complete 
parameter listing. The radial velocities and the derived solution 
are shown in Fig. 2. When fitting an eccentric radial-velocity 
curve the amplitude goes down to 56.8 km s~', for eccentricity 
e=0. 063(33). Throughout this paper we regard the orbit as cir- 
cular. 

Given the solution presented in Table [2] the orbital radius of 
the sdB star can be approximated by a s dB sin i = 11.5 Rq, which 
corresponds to a light-travel time of 53.6(1.6) seconds between 
orbital phases corresponding to closest and furthest distance to 
the Sun. 

The orbital solution combined with the mass function gives a 
lower limit for the mass of the companion of more than 0.63 Mq, 
if one assumes a canonical mass of 0.48 Mq for the sdB star. As 
the spectrum does not reveal clear evidence for light contribu- 
tion from a companion, it should be either an unseen compact 
object, or an unseen K star. As the 2MASS J = 15.38(5) and 
H= 15.35(9) magnitudes do not indicate a rising infrared flux 
that would reveal the presence of a K star in this system, we con- 
clude that the unseen companion is most likely a white dwarf. If 
the inclination of the system is smaller than i ^ 40 degrees then 
the companion may be a neutron star or black hole (z $ 25 de- 
grees). 
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Table 1. Log of the low-resolution spectroscopy of KIC 1 1558725 



Mid-exposure Date 


Barycentric JD 
-2455000 


S/N 


RV 

km s 


RV C1T 
km s ^ 


Telescope 


Observer/PI 


2010- 


-07- 


-28 00:24:52.0 


405.519273 


113.9 


-26.6 


8.6 


WHT 


JHT/CA 


2010- 


-07- 


-28 05:16:04.8 


405.721505 


94.6 


-30.6 


9.0 


WHT 


JHT/CA 


2010- 


-08- 


-12 10:05:19.3 


420.922376 


57.6 


-93.5 


10.0 


KP4m 


MDR,LF 


2010- 


-08- 


■12 10:16:03.9 


420.929837 


60.6 


-91.1 


9.8 


KP4m 


MDR,LF 


2010- 


-08- 


-13 03:41:29.7 


421.655828 


68.3 


-69.3 


7.4 


KP4m 


MDR,LF 


2010- 


-08- 


-13 03:52:04.9 


421.663178 


74.1 


-66.4 


8.1 


KP4m 


MDR,LF 


2010- 


-08- 


-14 08:24:41.8 


422.852489 


59.1 


-31.1 


8.0 


KP4m 


MDR,LF 


2010- 


-08- 


-14 08:35:42.3 


422.860134 


59.4 


-25.9 


8.5 


KP4m 


MDR,LF 


2010- 


-08- 


-14 08:46:01.9 


422.867305 


60.1 


-34.2 


6.9 


KP4m 


MDR,LF 


2010- 


-08- 


-25 22:54:52.3 


434.456682 


91.3 


-4.6 


2.1 


WHT 


RH0/CA 


2010- 


-08- 


-25 23:09:46.0 


434.467026 


106.2 


-17.4 


5.2 


WHT 


RH0/CA 


2010- 


-08- 


-27 00:36:26.9 


435.527210 


117.9 


-25.8 


9.2 


WHT 


RH0/CA 


2010- 


-08- 


-27 00:46:36.6 


435.534266 


112.3 


-31.2 


7.3 


WHT 


RH0/CA 


2010- 


-08- 


■27 21:56:54.7 


436.416408 


93.6 


-39.0 


8.2 


WHT 


RH0/CA 


2010- 


-08- 


-27 22:08:29.4 


436.424449 


90.5 


-63.8 


4.4 


WHT 


RH0/CA 


2010- 


-08- 


■28 23:18:34.2 


437.473103 


74.7 


-83.9 


8.3 


WHT 


RH0/CA 


2010- 


-08- 


■28 23:28:43.8 


437.480158 


71.5 


-90.5 


6.1 


WHT 


RH0/CA 


2010- 


-08- 


■29 21:53:35.3 


438.414074 


105.3 


-111.8 


9.1 


WHT 


RH0/CA 


2010- 


-08- 


■29 22:03:44.8 


438.421129 


107.7 


-128.2 


10.9 


WHT 


RH0/CA 


2010- 


-08- 


■31 01:31:47.4 


439.565587 


119.6 


-116.6 


6.0 


WHT 


RH0/CA 


2010- 


-08- 


■31 01:41:56.9 


439.572642 


109.0 


-117.0 


6.4 


WHT 


RH0/CA 


2011- 


-05- 


■31 01:56:05.9 


712.581494 


49.4 


-85.7 


18.4 


NOT 


JHT 


2011- 


-05- 


■31 04:48:22.9 


712.701140 


79.6 


-88.5 


14.7 


NOT 


JHT 


2011- 


06 


■01 02:35:14.4 


713.608708 


57.0 


-41.4 


17.1 


NOT 


JHT 


201 1- 


-06- 


■07 03-46-43 9 


719.658535 


67.2 


-123.2 


9.2 


NOT 


JHT 


2011- 


-06- 


■09 03:34:59.9 


721.650444 


49.6 


-99.0 


14.2 


NOT 


JHT 


2011- 


06 


■10 04:58:52.1 


722.708717 


43.9 


-103.3 


10.0 


NOT 


JHT 


2011- 


-06- 


■20 00:20:10.7 


/ Ji.J IJ'rJO 


1 


— / O. J 


1 A Q 


MOT 


TT4T 


2011- 


-06- 


■21 01:27:07.4 


733.561953 


36.3 


-69.3 


14.7 


NOT 


JHT 


2011- 


-06- 


■27 21:53:09.0 


740.413514 


56.4 


-148.6 


19.1 


NOT 


JHT 


2011- 


-07- 


■22 23:30:41.0 


765.481614 


38.1 


-22.1 


14.4 


NOT 


JHT 


2011- 


-07- 


■23 22:43:51.2 


766.449100 


49.2 


-28.4 


16.1 


NOT 


JHT 


2011- 


-08- 


■29 00:17:40.9 


802.514159 


42.1 


-110.0 


16.6 


NOT 


RO 


2011- 


-08- 


■30 00:38:59.5 


803.528944 


48.1 


-81.2 


16.6 


NOT 


RO 


2011- 


-08- 


■31 00:25:05.0 


804.519272 


38.2 


-50.1 


13.9 


NOT 


RO 



Assuming a radius for the sdB star of 0.2 Rq, and assuming 
a white dwarf companion, one may expect to see eclipses only 
if the inclination angle is higher than i ^ 88 degrees. We did not 
detect eclipses (see Fig. 3). 

2.2. Atmospheric parameters 

The spectra were shifted to remove the orbital motion, 
before being co-added to obtain high-S/N spectra (KP4m 
S/N =175; NOT S/N = 145; WHT S/N = 305) with minimal 
orbital line broadening, for all three observatories. We de- 
rive the atmospheric parameters of the star from each of 
these mean spectra, and produce a weighted mean using 

Table 2. Orbital solution of KIC 1 1558725 



system velocity [km s '] -66.1(1.4) 

radial-velocity amplitude K [kms -1 ] 58.1 (1.7) 

period P [day] 10.0545 (0.0048) 

phase [BJD - 2455000] 421.682 (0.062) 

X 2 34.48 

reduced^ 2 1.11 



the formal fitting errors as variance weights, as listed in 
Table |3] Our final adopted values are T eS = 2791 0(2 10)K, 
log g = 5 .4 1 0( 1 5) dex and log (N H JN H ) = -3.116(18) dex, which 
are quite compatible with the parameters T s s = 27400 K, 
log g = 5.37, log (Nhs/Nh) = ~2.8 found from the initial survey 
spectrum in lOstensen et al.l (2011). The errors on the adopted 
values are the errors on the weighted means, and reflect the 
spread of the individual measurements rather than the formal er- 
rors. Systematic errors related to model physics typically are of 
the order 500 K, 0.05, and 0.05, for 7/ eff , log g, and log (NnJN H ), 
respectively. 



Table 3. Atmospheric parameters of the subdwarf in 
KIC 11558725 



Telescope 






log(AWAW 




K 


cm s~ 2 




KP4m 


27950(110) 


5.378(20) 


-3.078(42) 


NOT 


27410(110) 


5.418(20) 


-3.124(32) 


WHT 


28010(50) 


5.417(11) 


-3.125(26) 



adopted 27910(210) 5.410(15) -3.116(18) 
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Besides the lines of the hydrogen Balmer series and Hei 
lines, clear lines of magnesium (448 1 A), and silicon (4552A, 
4567 A) are present. The high S/N mean WHT spectrum seen in 
the upper panel of Fig. Q]also allows the identification of a num- 
ber of lines from O n, N n as well as the C n doublet at 4267 A and 
the C in lines at 4647 and 4650A. By comparing the strength of 
these lines relative to LTE model spectra with explicit metal lines 
computed for the adopted parameters listed in Table [3] varying 
only the metallicity relative to solar, we find that oxygen is de- 
pleted to 4% of the solar abundance, nitrogen to 13%, and car- 
bon to 1.6%. Such high depletion is normal in sdB stars due to 
gravitational settling, and large deviati ons from the solar mix for 
individual elements are also common dHeber et al.ll2000l) . 



3. The orbital light curve from Kepler photometry 

We analysed the Kepler light curve of KIC 1 1558725 as obtained 
in quarters Q6-Q10, totalling roughly 15 mont hs of data with 
58.8 s sampling time ( i.e. short-cadence data; see Gillilan d et al.l 
120101) . The photometric time series spans BJD 2455372 to 
2455833, which is roughly the same range as our spectroscopic 
dataset. We analyse the light curve as transformed to fractional 
intensities AI/I. 

At first we used the raw light curves pro duced by the 
Keple r data processing pipeline (SAP_FLUX; iJenkins et alj 
2010). These light curves are extracted using a standard small 
pixel mask selected to minimise the noise in the output data. 
However, during the observations slight drifts were induced by 
guiding corrections and thermal effects in the focal plane, which 
implies that some light drifts in or out of the small pixel mask, 
causing trends that must be removed. For pulsations with much 
shorter timescales than the trends, this is not a problem, but with 
a 10 d orbital period in KIC 1 1558725, our concern was that the 
detrending could suppress the orbital beaming effect. We there- 
fore analysed also the raw Kepler pixel data using a large cus- 
tom pixel mask including all 20 to 23 pixels with target signal, 
as opposed to 4 to 6 pixels in the small pixel mask of the stan- 
dard pipeline extraction. This recovers 36% more flux from the 



target and the monthly chunks can be detrended with simple lin- 
ear trends. This approach gave a beaming amplitude that is 34% 
higher than our first estimate from the pipeline-reduced data. 
However, these large-aperture light curves are substantially more 
noisy and contain many more bad points. Hence we only use the 
raw, custom-mask pixel data to determine the amplitude of the 
Doppler beaming, but use the standard-mask pipeline-reduced 
raw data for our pulsation-frequency analysis in Sect. 4. 

A Fourier analysis of the Kepler light curve reveals the or- 
bital period to be 10.0516(27)d, in perfect agreement with the 
period found from spectroscopy. In Fig. [3] we show the Kepler 
light curve folded on the orbital period into 50 phase bins. A 
sine fit to the folded light curve, leaving only the amplitude as a 
free parameter, gives an amplitude of the orbital effects seen by 
Kepler of 238(6) ppm, and this amplitude is consistent with that 
derived from the Fourier analysis. The fit is shown as a solid line 
in the figure. We do not find orbital harmonics. 

3.1. Doppler beaming 

The high precision Kepler data permits us to accurately explore 
the low-amplitude Doppler beaming effect, something that is 
very hard to do with ground based data. This effect is induced 
by stel lar motion in a binary orbit and causes brightness modu- 
lation (Rvbicki & Lightman 1979). The Doppler beaming effect 
permits an estimate of the radial velocity without resorting to 
spectroscopic data. This effect was d e tected for the first time 
in Kepler data by Ivan Kerkwiik et all (1201 Oh and fo r a plane- 
tary system in Caftordatabv lMazeh & Faiglerl OOlO). A confir- 
mation of the correspondence between Doppler beaming ampli- 
tudes and radial velocities was first established bv lBloemen et al.l 
(1201 lh . 

The light curve of KIC 11558725 displays a brightening 
of the sdB star at the orbital phases where the star is ap- 
proaching us in its orbit (i.e. when its orbital radial veloci- 
ties are negative, see Fig. [3), and this is exactly the effect ex- 
pected by Doppler be aming. Note that unlike KPD 1946+4330 
ilBloemen et1d]|201 lb . KIC 11558725 does not show any sign 
of ellipsoidal deformation which in the closest sdB+WD bina- 
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ries produces a str ong harmonic signal at P rb/2 (see also e.g. 
ISilvotti et al.ll2012h . This is consistent with the much longer or- 
bital period of KIC 11558725 as opposed to the 0.4 d period of 
KPD 1946+4330. 

In the case of Doppler beaming the observed flux from the 
target, F^, is related to the emitted spectrum and orbital velocity 
as 

Fi = F %l [l-B V -^ (1) 

and the beaming amplitude relates to the orbital radial-velocity 
amplitude K as 

Ab D = B— , (2) 

c 

where As is the amplitude of the beaming signal in the light 
curve (see Fig. 3), and D the decontamination factor discussed 
below. 

The beaming factor B contains three terms. For an approach- 
ing source a contribution of + 1 comes from the enhanced photon 
arrival rate, and another term of +2 from the geometrical aber- 
ration of the wavefronts. The last term comes from the Doppler 
shift of the spectrum, which for a hot sdB star produces a neg- 
ative effect as the blue-dominated flux gets blue-shifted out of 
the observed Kepler bandpass when the star is app roaching. We 
comp ute this last term using TMAP model spectra (I Werner et alj 
l2003h that cover th e entire Kegler bandp ass, following the pro- 
cedure described in Bloemen et 

With the beaming factor B- 1.403(5), and the spectroscop- 
ically derived orbital radial-velocity amplitude, we compute a 
predicted total beaming amplitude for the Kepler bandpass of 
270ppm. The observed beaming amplitude of Afi=238ppm is 
88% of the predicted value. 

The Kepler pixels onto which our target is imaged suffer 
from contamination from neighbouring objects, and from pass- 
ing charge from brighter sources when clocking out the CCD. 
According to the Kepler Target catalogu^l the contamination 
factors for KIC 1 1558725 are 0.098 (Q6,Q10), 0.1 10 (Q7), 0.079 
(Q8) and 0.106 (Q9), or 0.0982 on average, implying that all 
periodic amplitudes derived from these Kepler fluxes should be 
multiplied by D- 1/(1-0.0982)= 1.109 to get the intrinsic am- 
plitudes of KIC 1 1558725. When applying this decontamination 
factor we find that the photometric amplitude of the Doppler 
beaming as seen by Kepler is consistent with the spectroscop- 
ically derived radial-velocity amplitude within the errors of the 
data. Given the observed amplitude Ag, we use Eq. (2) to de- 
rive a value of ,K=56.5(1.4)kms~ 1 for the orbital radial-velocity 
amplitude, from Doppler beaming. 

The fact that this value is consistent with that of the spec- 
troscopic value proves that the companion of the sdB does not 
significantly contribute to the observed Doppler beaming, con- 
sistent with a compact nature of the companion. 

4. The pulsation spectrum 

We are interested in the pulsation frequencies in KIC 1 1558725 
for three reasons. Firstly, we intend to use the pulsations as 
clocks to derive the orbital light-travel time delay and conse- 
quently another independent determination of the orbital radial- 
velocity amplitude. Secondly, the pulsations may reveal the in- 
ternal structure of the subdwarf through a detailed seismic study, 
for which we intend to derive mode identifications in Sect. 6. 

1 http://archive.stsci.edu/kepler/kepler_fov/search.php 



Furthermore, the frequency splittings may disclose the rotation 
period of the subdwarf. 

To investigate the pulsational frequency spectrum we used 
the Kepler light curves as produced by the standard Kepler 
data processing pipeline. The fractional object intensities AI/I of 
quarters Q6-Q10 were first detrended, and subsequently outly- 
ing points were discarded leaving 625617 data points with stan- 
dard deviation cr=2600ppm . An iterative prewhitening process, 
involving a standard fast-fourier transform to find peaks and sub- 
sequent non-linear least-squares (NLLS) sine-curve fits to sub- 
tract the pulsational content mode by mode from the original 
data, was used to derive the frequency list of Table 4. In total we 
extracted 166 pulsational frequencies, in addition to the orbital 
frequency. The standard deviation of the data after removal of all 
significant periodic signals was <x= 1900ppm. 

The standard deviation in the Fourier amplitude spectrum 
of the original data amounts to cr FX =4.2ppm, and we adopt 
4cr FT =17ppm as the threshold of significance of the peaks in 
the Fourier amplitude spectrum. 

We note that the observed amplitudes in Table 4 have not 
been corrected for the Kepler decontamination factor D (see 
above), nor for the amplitude smearing due to the effective expo- 
sure time of 58.8 seconds. Both effects cause the observed am- 
plitudes to be smaller than the intrinsic amplitudes, and the latter 
affects mostly the shortest periods (p-modes). For the strongest 
pulsation, with a period of 3641.1 s, the exposure time leads to 
a decrease of the amplitude of 4.3%, while for the shortest ex- 
tracted period of 197.9 s the amplitude decrease is 14%. 

To illustrate the fact that all pulsational frequencies are 
present throughout the Kepler run, we show a section of the 
Fourier transform in a dynamic form in Fig. 4. There is clear 
beating among the stronger frequencies; in fact, these beatings 
show up as resolved frequencies in the Fourier transform of the 
full Kepler data set, and may be attributed to the rotation of the 
subdwarf (Sect. 7). 

4.1. Combination frequencies 

We do not find evidence for harmonics of the pulsation frequen- 
cies. We do, however, find evidence for combination frequen- 
cies. When computing the residual frequency Sf = fa — fi — f\ 
for all combinations of 168 frequencies in Table 4, we find 13 
combinations that are consistent with a residual 6f-Q within the 
errors of the frequencies. Thirty-five different frequencies are in- 
volved in these 13 combinations. One of the combinations is an 
orbital alias discussed in Sect. 5.1. See Table 5 for the list of 
frequency combinations. 

To investigate the significance of this, we take our observed 
frequencies and perturb each of them by a random offset which 
is big enough to upset real combinations, but small enough such 
that the distribution of these randomised frequencies looks like 
that of the observed distribution. We perturb the frequencies by 
+ 1 /vHz. For each randomised frequency we adopt the frequency 
error of the original unperturbed peak. For 50 sets of such ran- 
domised frequencies, we find an average of 3.7 combinations 
that satisfy 6f-Q within the errors of the frequencies. 

Hence, for the set of 13 combinations that we find in our 
168 frequencies, we expect that about 4 are by-chance com- 
binations. We assume that the combinations with the lowest- 
amplitude peaks are most likely to be by-chance combinations, 
as these peaks have the largest error on their frequency. These 
probable frequency combinations are marked in Table 4, and we 
take care when assigning these modes in period spacings and 
frequency multiplets. 
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Fig. 4. A dynamic Fourier spectrum from quarters Q6 through Q10 computed for 20 day stretches of data. The pronounced beatings 
show up as resolved frequencies in the Fourier spectrum of the complete data set. The horizontal lines reflect data gaps. 



5. Using the pulsations as clocks to derive the orbit 

In analogy to the effect that the Earth's orbit has on arrival times 
of a variational signal, and on the apparent frequencies of that 
signal, the orbit of the sdB alters the phases and frequencies of 
its pulsations as perceived by a distant observer. 

As we already established in Sect. 2, the light-travel time 
between conjunctions is 53.6 s, which amounts to a significant 
phase change for the 300 second p-mode pulsations in this 
star. Similarly, the Doppler frequency change of the pulsations 
at orbital quadratures amounts to Af = f * K/c, which, with 
/f=58 kms -1 , even for g-modes gives frequency shifts of 0.02- 
0.2 yuHz. These are certainly measurable with our 460 day data 
set with frequency resolution of 0.025//Hz. If the orbit is known, 
one can simply correct for these effects by converting all obser- 
vation timings to the centre-of-mass frame of the system, equiv- 
alent to the common approach of conversion of observation tim- 
ings to the barycentric frame of the solar system. 

If th e orbit is not known, o ne can use the pulsations as clocks 
(see e.g. lHulse & Tavlorll975b to derive the light-travel time that 
corresponds to the radius of the orbit, which relates to the radial- 
velocity amplitude K that can be obtained from spectroscopy or 
the Doppler-beaming curve as 



flsdB sin i 



(3) 



c c 2n 

where we introduce the R0mer delay, A?r, to represent the light- 
travel time. 

For a circular orbit, the light-travel delay as a function of the 
subdwarf's position in its orbit can be written as 



7\lek.><'> = A/r COS I ^-(*- Toft) 



(4) 



where T Q1 ^ is the time at which the subdwarf is closest to the Sun 
in its orbit, corresponding to the orbital phase listed in Table 2. 

We did not detect any second- or third-harmonic peaks, nei- 
ther for the orbit nor for the pulsations. Hence, for sinusoidal 



signals the Kepler light curve of KIC 1 1558725 can be approxi- 
mated by 



A/(f) / In 

—— = Ab sin I — — (t — Tixb) 

\ *orb 

+ Ai 'P uls Sin ( 27r ^i,puls(f - T^puls + r delay (f))) , (5) 



where the first term describes the orbital beaming effect, and 
where all individual pulsations are affected by the same orbital 
light-travel delay T^ e i ay (t), Here, the phase of the individual pul- 
sations, ri,p U i s , is expressed in the time domain rather than as an 
angle. Note that the above sum of sine curves is equivalent to the 
model that we fit as part of the prewhitening procedure described 
in Sect. 4, with the addition of a phase delay that introduces just 
one extra parameter, i.e. the amplitude AfR of the light-travel de- 
lay. 

We apply the above Eq. (5) to derive a third independent 
measurement of the radial-velocity amplitude of the subdwarf in 
KIC 11558725. For this purpose, we included the above phase 
delay in our NLLS Levenberg-Marquardt routine, and fitted the 
above model (5) to the Kepler Q6-Q10 light curve. In the model 
we included the 70 strongest frequencies that are identified to 
be part of t-\ and 1-2 sequences of the subdwarf (see next sec- 
tion), as to avoid possible pulsations from the unseen component 
in the binary. 

Simultaneously fitting all 70 pulsational amplitudes, fre- 
quencies and phases, together with the R0mer delay and 
phase as free parameters, we find A?r= 26.5(1.5)s and r or b= 
421.744(88) [BJD - 2455000], which is in perfect agreement 
with the spectroscopic parameters in Table 2, and the Doppler 
beaming amplitude. Given the observed R0mer delay A?r, we 
derive a value of /T=57.5(3.2)kms _1 for the orbital radial- 
velocity amplitude, from the light-travel time effect. 
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Fig. 5. a: The pulsation spectrum KIC 1 1558725 as observed (upper panel), and after correcting for the R0mer delay (lower panel) . 
The line in the upper panel indicates our 4-cr detection limit, b: The amplitude of five pulsation peaks as a function of R0mer 
correction, c: Same as a), but for a model light curve with the same duration and sampling frequency as the data in a), d: The 
predicted amplitude of the central peak and first alias as computed for a pure sine function distorted by an increasing R0mer delay. 
The upper curve indicates the total power of the central peak and the aliases. This sum departs from unity when the unaccounted-for 
second orbital alias pair starts to contribute. Entries in the legends have the same top-to-bottom order as the curves. 



5.1. Orbital aliases in the Fourier domain 

Above we have solved the orbit by directly fitting the light curve 
with a parameterised model of the pulsations that includes the or- 
bital light- travel delay. A differ e nt app roach was recently inves- 
tigated by IShibahashi & Kurta (I2012I) . who demonstrated how 
one can derive the orbit from the ratio of the amplitude of the 
central peak to the orbital aliases that the R0mer delay intro- 
duces in the observed Fourier spectrum. 

In the case of KIC 11558725 we find that, of our 166 
extracted pulsation frequencies, only a few of the strongest 
g-modes and a few of the highest-frequency pulsation modes in 
the weak p-mode regime have orbital aliases strong enough to 
be close to the amplitude level of significance. 

Around the strongest g-modes the orbital sidelobes appear 
in frequency regions populated by weak unresolved structure 
and prewhitening residuals. A more clear example of the or- 
bital sidelobes is given by the two peaks at 3073.5 and 3075.0 
//Hz, shown in the top panel of Fig. 5a. The sidelobes of each 
of the two central peaks can just be discerned close to the 4.0-cr 
line at a separation of F or \, = 1.15 //Hz. Since the amplitudes of 
these alias es are hardly significant, di rect application of the ratio 
method of Shibah ashi & Kurtzl (120121) leads to inaccurate results 
for this particular dataset. 

The two central peaks, however, are at ~16 <x, and their am- 
plitudes can be measured with very high precision. We can there- 



fore exploit the energy theorem in Fourier analysis that implies 
that the total power in the Fourier transform should not change 
due to the time-distorting orbital motion. All the power observed 
in the sidelobes will therefore be restored into the central peak 
if the light curve is corrected for the light-travel time. In the bot- 
tom half of Fig. [5^ we show that this is indeed the case. When we 
perform a light-travel time correction on the observation timings 
of the light curve, in the same way as the barycentric correction 
is done to remove the spacecraft orbit, then the orbital aliases in 
the Fourier domain disappear and the power is transferred to the 
central peak. For comparison, Fig. |5J; shows a simulation using 
noise-free, evenly-sampled data modelled with and without an 
imposed R0mer delay. 

We can also use the amplitude of the central peak to measure 
the R0mer delay, and get another estimate of the radial-velocity 
amplitude. If we take our observed light curve and adjust the ob- 
servation timings for various values of the R0mer delay, we find 
that the Fourier amplitudes of the pulsation peaks vary as shown 
in Fig.[5t>. Clearly, all five modes reach maximum amplitude for 
roughly the same delay correction. The delays and peak ampli- 
tudes are given in Table [6] and if we take an average using the 
amplitudes as weights we find Af« = 26.4(0.7) s. The gain in am- 
plitude for these modes is about 6%, when the optimal R0mer 
delay is applied. 
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A model of how the observed amplitude of a pulsation peak 
taken to reside at 3073.5 pMz drops with increasing Af# is shown 
in Fig.[5Jl. The corrected amplitude of the central peak is found 
to be 6.1% higher in a simulated light curve with no R0mer de- 
lay than in one with a 26 s delay, as observed. The predicted 
amplitude of the orbital alias is found to be 25.2%. Rearranging 
Eq. (27) in ShibahashL&Kurt^ J2012h in terms of this ratio, we 
find 

A +1 +A_! nAt R 

74 = ~P~ ' (6) 

^-^O r osc 

which corresponds to 24.2% for Af# = 25 s and f osc = 325 s, 
which demonstrates that our numerical modelling (Fig. 5d) is 
consistent with their analytical derivation. 

In Fig. 5d we show that the total power contained in the cen- 
tral peak and its aliases is preserved while applying observation- 
timing corrections. The curve drops off for large delays as we 
have only accounted for the first pair of a full comb of orbital 
aliases. 

6. Pulsational period spacings 

Recently, Reed et al] d201 lb have revealed that g-modes in sdB 
stars show sequences with almost fixed period spacings. This is 
contrary to what models predicted up to now, i.e. strongly vari- 
able period spacings caused by mode trapping due to a highly 
stratified internal structure. The observed regular period spac- 



ings indicate that internal mixing processes must be consider- 
ably stronger than presumed. As the period spacings follow the 
asymptotic relation for g-modes, they change with the value of 
the spherical-harmonic degree I of the modes. Hence, the spac- 
ing s allow us to deter mine the /"-value of the modes directly. 

iReed et al.l (1201 ll) already identified t values for the pul- 
sations that could be resolved in the Kepler survey data of 
KIC 11558725. They identified an l=\ sequence with period 
spacing AIi£ = i =246.8 s, and an 1-2 sequence with period spac- 
ing AIT = 2=142.6 s. 

Here we use the same methods of IReed et alj d201 lb to de- 
rive the period spacings, using the much longer Q6-Q10 data 
set. First we defined the period spacing for l=\, as marked in 
Table 4. This was achieved iteratively. As 1=1 modes suffer 
the least geometric cancellation (see e.g. lReed et al.l l2005). and 
hence should have highest amplitudes if all modes are driven at 
a similar level, we started off with the highest-amplitude modes. 
Subsequently we worked with an €=\ sequence consisting only 
of doublets, and finally all the l=\ periods marked in Table 4. 
For each step we derived mutually consistent values o f the period 
spacin g. We recover all the £-\ identifications from Reed et al. 
(1201 ll) except for two periods not detected in the Q6-Q10 data. 

Including all 57 l-\ multiplet candidate peaks from 
Table 4 in a linear regression P(n) = C + n * Anf = i, we find 
An f= i=248.87(18)s with RMS=19s. Here, the mode sequence 
number n should be equal to the radial number of nodes in the 
pulsation, n', but shifted by a fixed offset that we cannot deter- 
mine. When fitting for each of the 29 multiplets (mostly dou- 
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blets) only the peak closest to this first fit, we obtain a consistent 
period spacing AII^i =248.68(23) s, with RMS=14s. 

Given that for the asymptotic g-mode domain we expect 
AIl{ = 2=AIl{ = i/ V3, we used a spacing of AIT = 2=143s as a 
starting point for the 1-2 sequence. Again we iterated adding 
the highest-amplitude peaks and some frequency multiplets that 
cannot be l-\ modes. Including all 76 (=2 multiplet candidate 
peaks in a linear regression, we find AIT = 2=143.44(08)s with 
RMS=16 s. When fitting for each of the 34 multiplets only the 
peak closest to this initial fit, we obtain a consistent period spac- 
ing AI^ =2 =143.37(10)s, with RMS=14s. 

Of 166 pulsation frequencies, 23 have periods shorter than 
900 s which we did not include as g-mode period-spacing can- 
didates. From the 143 g-mode pulsations we find that 108 can 
be matched with either the t=\ or 1-2 sequence, when accept- 
ing periods deviating up to ~35 s from the linear regression and 
allowing for any possible value of the azimuthal order m. We 
plot echelle diagrams for the (=1 and 1-2 sequences in Fig. 6, 
and list all t=\ and £-2 matches in Table 4. Periods below 900 
seconds that do match any of the 1= 1 or 1-2 sequences are indi- 
cated as such in Table 4 for completeness. 

Obviously, the period spacing in KIC 11558725 is not very 
strict, which makes conclusive identification of individual modes 
difficult, as by-chance identifications of modes that do not be- 
long to a sequence are probable. However, the fact that the spac- 
ings vary within a sequence may hold information about the in- 
ternal structure that can be derived from detailed seismic mod- 
elling. From Fig. 6 we find that between sequence numbers 8 
and 18 in the t=\ sequence, corresponding to pulsation periods 
between 3300-5900 s, there is a remarkable bump away from 
the average spacing. Whereas the t=\ sequence shows coher- 
ence over a large period range, the 1-2 sequence is dominated 
by a clumping of multiplet candidates in the period range 2900- 
3700 s. 

Pulsation modes that could erroneously be accepted in any of 
the sequences may stem from any of the following: modes with 
higher degree t that are strong enough to be picked-up by Kepler, 
modes in combination frequencies; modes that are tidally forced 
to have large amplitudes; modes that originate in the companion. 
If a mode in a combination frequency is due to non-linear or 
resonance effects between two other modes, its assignment to 
a period-spacing sequence may not be justified. Therefore we 
mark the modes that take part in frequency combinations (see 
Table 5) with a circle in Fig. 6. 

Tidally forced oscillations can have large amplitudes when 
close to an eigenfrequency of the star, and can in principle be of 
any degree (, although an 1-2 excitation can be considered the 
most probable, in which case tidally induced pulsations should 
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Table 6. Measured R0mer delay from high-frequency pulsa- 
tions. 



Frequency 


At R 


A 


iMz 


s 


ppm 


3059.896 


26.04(05) 


62.2(1) 


3073.467 


25.45(05) 


66.2(1) 


3075.021 


26.25(07) 


53.6(1) 


3102.271 


29.35(15) 


31.4(2) 


3103.849 


29.23(13) 


29.9(1) 



Notes. The amplitude A is the recovered amplitude when applying the 
corresponding R0mer delay At R . 



Fig. 7. Selected multiplets. The bottom 5 panels show the 
~0.13juHz splitting seen in some of the t-\ doublets. The 
strongest pulsations we detect in KIC 1 1558725 are amongst the 
l- \ peaks shown here. The top panels show sample regions with 
sequences of frequencies with splittings in the 0.37-0.58//Hz 
range. Amplitudes are given in ppm. The significance cut-off 
amplitude is 17 ppm. 

be part of the 1-2 sequence. The error on the orbital frequency is 
not small enough to exclude some of the strongest pulsations in 
Table 4 from being a multiple of the orbital frequency. We note 
however, that for the observed frequencies in KIC 11558725 
tidal forcing would occur at very high harmonics of the orbital 
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0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 0.6 
Frequency splitting [u.Hz] 

Fig. 8. Histogram of splittings listed in Table 4. The blue boxes 
reflect splittings from multiplets (mostly doublets) that were 
matched with the l=\ sequence in Sect. 6. 

frequency (i.e. >300 for the strongest pulsation frequencies), 
which may diminish the forcing efficiency. Nevertheless, a re- 
cent case of tidally driven pul sations at high orb ital harmonics 
(90th and 91st) is presented bv lWelsh et aiT(l201 lh . 

Even though pulsating white dwarfs are rare as the instability 
strips are small, possible pulsations in the supposed white-dwarf 
companion could be picked up by Kepler. Assuming a white 
dwarf with a large radius 0.013 Rq (log(g)=8.0, M=0.63 Mq), 
with a similar or somewhat cooler temperature as the subdwarf, 
and with intrinsic pulsation amplitudes at the 10-20% level, we 
conclude that such a white dwarf would contribute pulsation fre- 
quencies at observable amplitudes to the frequency spectrum of 
KIC 11558725, even if the white dwarf is not detected in our 
spectra, or in the available broad-band imaging. Hence we can 
not exclude this possibility, when trying to match period spac- 
ings or frequency splittings in KIC 1 1558725. 

If any of the significant peaks we have detected were to orig- 
inate from the white dwarf, their amplitudes would not be max- 
imised by shifting to the rest frame of the sdB, as we demon- 
strated is the case for the short period peaks in Sect. 5.1. Rather, 
a R0mer delay with the opposite sign and an amplitude corre- 
sponding to the radial velocity of the white dwarf would have 
to be imposed to recover the full amplitude. Using the tech- 
nique of Sect. 5.1 we found only one peak in Table 4 that shows 
an orbital phase behaviour that is opposite to that expected for 
the subdwarf, and hence could originate in the supposed white- 
dwarf companion: it is the peak with the shortest periodicity at 
5053.5 //Hz. As this peak is just at the detection limit we cannot 
rule out that it is just a random effect. If it persists in the next year 
of Kepler observations, it might turn out to be the first indication 
of an sdB+WD system where both stars are pulsating. 

7. Frequency splittings and the rotation period 

Fixed frequency splittings are expected as a consequence of ro- 
tation. To the first order in the rotation frequency £2, the observed 
frequencies of the modes are altered by mQ(l - C„>f), with m the 
azimuthal quantum number of the spherical harmonic. Hence, 
for non-radial modes of given degree I, we expect frequency 
multiplets of 2^+1 peaks to occur. However, not all peaks of a 
full multiplet need to be excited, and, furthermore, cancellation 



effects for a given inclination angle can favour the observability 
of som e modes over others within a multiplet (see e.g. lReed et al.l 
l200l . 

In the case of KIC 1 1558725 the frequencies mostly do not 
appear as full multiplets, which makes the interpreting of the 
observed splittings difficult. See Table 4 for the observed split- 
tings. We find many splittings in the range of 0.1-0.5//Hz. In 
Fig. 7 we show two sample frequency sequences that show a 
frequency splitting in the 0.4-0.55 //Hz range, that seems typical 
for frequency groups that do not match the {-I period spacing. 
We also show in Fig. 7 a series of doublets that match the l=\ 
period spacing, which all have a frequency splitting on the or- 
der of 0.1 3 //Hz, and which in fact cause the prominent beating 
seen in Fig. 4. We also detect some doublets with spacings ap- 
proximately twice this (~0.23 //Hz) and we interpret the smaller 
frequency spacings as consecutive m while the larger spacings 
would be the \m\-\ pairs. In Fig. 8 we present a histogram of 
the splittings listed in Table 4. 

We tentatively identify this splitting of 0.13 //Hz as the rota- 
tion frequency of the subdwarf in KIC 1 1558725. We regard this 
result as inconclusive as other doublets that seem to match the 
l-\ period spacing show a larger frequency splitting (>0.3 //Hz). 
We note however that some of the peaks in these i=\ doublets 
also match the 1-2 period spacing, and hence could be 1=1 
modes instead. 

Assuming that the Ledoux constant C n >i can be approximated 
by l/((((+lJ) as in the asymptotic g-mode regime, the observed 
frequency splitting of 0.1 3 //Hz in the £=l doublets leads to a 
rotation period of the subdwarf of 45 d. This consequently im- 
plies that the rotation of the subdwarf is not phase-locked with 
the orbit. 

Even if we assume that the larger frequency splittings, in 
the 0.4-0.55 //Hz range, belong to multiplets of i> 1 modes and 
are due to rotation, we would still conclude with certainty that 
the subdwarf rotates subsynchronously with respect to the or- 
bit. Here the case of a splitting of 0.55 //Hz for a supposed 1=2 
multi plet provides the lower limit to the rotation period of 18 d. 

iPablo et ail d201 ll 1201 2|) discuss the frequency splittings in 
the three of the four known Kepler sdB+dM systems, which all 
three have orbital periods on the order of ten hours. From the 
frequency splittings they derive that for each of their cases the 
subdwarf rotation is subsynchron ous w i th resp ect to the orbit. 
Following their description of the IZahnl (1 19751) synchronization 
timescale, we expect for KIC 1 1558725 a much longer synchro- 
nization timescale than for those sdB+dM systems, as the orbital 
period is much longer and the synchronization timescale scales 

17 /3 

with PJ . In fact, we expect that the rotation period we measure 
now is a good indicator of the rotation period of the subdwarf af- 
ter it just settled on the extreme horizontal branch. 

8. Summary and Conclusions 

From our new low-resolution spectroscopy we discovered that 
KIC 11558725 is a binary consisting of a B subdwarf and an 
unseen companion, likely a white dwarf. We found a circu- 
lar orbit with P or b= 10.05 d, and a radial-velocity amplitude 
of 58km s . From the high signal-to-noise average spectrum 
we redetermined the atmospheric parameters of the subdwarf: 
r eff = 27910K and log g = 5 .4 1 dex. 

Five quarters of short-cadence Kepler data of KIC 1 1558725 
reveal Doppler beaming at the 238 ppm level. The amplitude of 
the Doppler-beaming modulation in the light curve corresponds 
to a radial-velocity amplitude that is fully consistent with that 
determined from the spectroscopic data. 
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We developed a new method that uses the pulsations as 
clocks to measure the orbital light-travel effect, or R0mer de- 
lay, directly from the Kepler light curve. The measured delay 
is again in perfect agreement with the spectroscopically deter- 
mined radial-velocity amplitude. This method allows the light- 
travel delay to be measured on many pulsations simultaneously, 
allowing the determination to be accurate, and hence can be use- 
ful for stars that have many weak pulsation frequencies such as 
sdB stars. 

For the 300 second /?-modes in KIC 1 1558725 the pulsation 
amplitude varies slightly when transforming the observation tim- 
ings to that of the center-of-mass of the binary. This effect will 
only show up for short-period pulsations in a long orbit. We have 
been able to measure the R0mer delay from this effect. 

One may wonder if it is still necessary to use spectroscopy 
to solve orbits for targets in the Kepler field. For the case of 
KIC 11558725 the phasing of the spectroscopic results with re- 
spect to the Kepler light curve proved essential in order to estab- 
lish that the light-curve modulation is due to Doppler beaming, 
and not due to reflected light from the companion or to a contam- 
inating object in the course Kepler pixel aperture. In general, the 
Doppler beaming amplitude may not reflect the radial-velocity 
amplitude in case the companion significantly contributes to the 
observed combined beaming amplitude. However, the phase and 
amplitude of the light-travel delay as determined by using the 
pulsations as clocks provide the same orbital constraints as de- 
rived from spectroscopy. 

We extracted 166 pulsation frequencies from the Kepler light 
curve, among which 13 are found to be possible combination 
frequencies. We found 6 /?-modes, and 160 frequencies in the g- 
mode domain, demonstrating the potential for a seismic analysis 
of this star. 

We showed that many of the pulsation frequencies match pe- 
riod spacings of (=1 and 1=2 sequences, implying that the in- 
terior structure of the subdwarf cannot be heavily stratified. The 
period-spacing sequences will aid in the identification process of 
the modes in a future seismic study of this object. 

Although we detect many pulsation frequencies, we see lit- 
tle evidence for complete multiplets. The frequency splittings 
are many in KIC 1 1558725 and are difficult to fit in a standard 
picture. We attribute the smallest splittings that we find to l=\ 
doublets, and derive a rotation period of 45 d. As a lower limit 
for the rotation period, we derive a value of 18 d, which implies 
that the subdwarf rotates subsynchronously with respect to the 
orbit. 

KIC 1 1558725 is the first sdB pulsator in the Kepler sample 
with a confirmed compact stellar-mass companion. The 10 day 
orbit is in the long end of the period range of the ~100 known 
sdB binaries that have periods compatible with the common- 
envelope ejection scenario. Assuming a canonical sdB mass of 
0.48 Mq, we derived a lower limit for the mass of the compan- 
ion of 0.63 Mq. The distance between the two companions is 
^20 Rq, which implies that if the sdB is a result of a common- 
envelope phase the progenitor must have been close to the max- 
imum radius for a red giant, near the tip of the red-giant branch. 
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Table 4. Extracted pulsational frequencies. 



Frequency 


Period 


Amplitude 


S/N Splitting n (=l 


n l=2 




s 


ppm 


//Hz 





5053.4659 


(20) 


197.88399 


(0.00008) 


23 (4) 


5.4 










3103.8490 


(17) 


322.18062 


(0.00017) 


28 (4) 


6.5 


1.58 








3102.2724 


(16) 


322.34436 


(0.00016) 


29 (4) 


6.8 










3075.0214 


(09) 


325.20099 


(0.00010) 


50(4) 


11.7 


1.55 








3073.4671 


(07) 


325.36545 


(0.00008) 


63(4) 


14.8 


1.15 








3072.3162 


(24) 


325.48733 


(0.00026) 


19(4) 


4.4 








OA 


3059.8960 


(08) 


326.80849 


(0.00008) 


59(4) 


13.8 










1463.3977 


(25) 


683.34123 


(0.00116) 


18(4) 


4.2 






-1 




1444.0731 


(19) 


692.48572 


(0.00091) 


24(4) 


5.6 








* 


1422.9117 


(11) 


702.78430 


(0.00053) 


43 (4) 


10.1 










1420.4369 


(17) 


704.00875 


(0.00083) 


27 (4) 


6.3 


0.24 








1420.1919 


(22) 


704.13022 


(0.00111) 


20 (4) 


4.7 








* 


1399.2750 


(23) 


714.65579 


(0.00119) 


20 (4) 


4.7 










1396.5427 


(21) 


716.05400 


(0.00109) 


22(4) 


5.1 


0.13 








1396.4129 


(20) 


716.12059 


(0.00105) 


23 (4) 


5.4 


0.27 








1396.1397 


(20) 


716.26072 


(0.00104) 


23 (4) 


5.4 








* 


1370.4145 


(11) 


729.70621 


(0.00057) 


43 (4) 


10.1 








* 


1308.1017 


(23) 


764.46658 


(0.00136) 


20(4) 


4.7 








* 


1292.1409 


(13) 


773.90943 


(0.00077) 


35 (4) 


8.2 








* 


1270.2485 


(13) 


787.24755 


(0.00083) 


34 (4) 


8.0 










1265.4160 


(19) 


790.25394 


(0.00116) 


25 (4) 


5.8 










1237.3675 


(20) 


808.16736 


(0.00134) 


22 (4) 


5.1 











1139.1458 


(23) 


877.85072 


(0.00175) 


20(4) 


4.7 










1128.9600 


(21) 


885.77098 


(0.00162) 


22 (4) 


5.1 










1083.6435 


(26) 


922.81273 


(0.00220) 


18(4) 


4.2 




-2 




* 


1042.2709 


(20) 


959.44349 


(0.00188) 


22 (4) 


5.1 


0.53 




T 




1041.7417 


(12) 


959.93086 


(0.00115) 


37(4) 


8.7 


0.26 








1041.4816 


(20) 


960.17062 


(0.00183) 


23(4) 


5.4 


0.27 




1 




1041.2151 


(14) 


960.41633 


(0.00131) 


32(4) 


7.5 


0.52 








1040.6927 


(11) 


960.89840 


(0.00104) 


41(4) 


9.6 


0.54 








1040.1546 


(20) 


961.39554 


(0.00187) 


23 (4) 


5.4 


0.58 








1039.5769 


(08) 


961.92981 


(0.00073) 


58(4) 


13.6 






1 




915.9922 


(20) 


1091.71231 


(0.00240) 


23 (4) 


5.4 






2 




888.0072 


(20) 


1126.11695 


(0.00248) 


23 (4) 


5.4 










875.0949 


(18) 


1142.73315 


(0.00237) 


25 (4) 


5.8 




-1 




* 


826.2266 


(13) 


1210.32167 


(0.00189) 


35 (4) 


8.2 










804.6305 


(09) 


1242.80646 


(0.00134) 


53(4) 


12.4 


0.44 




T 




804.1940 


(26) 


1243.48108 


(0.00408) 


17(4) 


4.0 


1.07 




3 




803.1235 


(14) 


1245.13844 


(0.00219) 


32(4) 


7.5 


0.56 








802.5678 


(17) 


1246.00067 


(0.00267) 


27 (4) 


6.3 






1 




774.4270 


(13) 


1291.27732 


(0.00217) 


35 (4) 


8.2 










741.5079 


(26) 


1348.60322 


(0.00468) 


18(4) 


4.2 










729.2202 


(24) 


1371.32792 


(0.00447) 


19(4) 


4.4 










712.5805 


(03) 


1403.35031 


(0.00067) 


135 (4) 


31.7 


0.40 


T 


T 


* 


712.1802 


(07) 


1404.13908 


(0.00137) 


66 (4) 


15.5 


0.41 








711.7661 


(05) 


1404.95599 


(0.00098) 


92(4) 


21.6 


0.37 




4 




711.3935 


(05) 


1405.69176 


(0.00104) 


87 (4) 


20.4 


0.34 









711.0567 


(09) 


1406.35756 


(0.00183) 


49 (4) 


11.5 




1 


i 




685.7967 


(20) 


1458.15803 


(0.00431) 


23 (4) 


5.4 










683.9097 


(23) 


1462.18144 


(0.00502) 


19(4) 


4.4 










663.5667 


(11) 


1507.00754 


(0.00253) 


41(4) 


9.6 


0.52 








663.0469 


(19) 


1508.18887 


(0.00439) 


24(4) 


5.6 










659.0236 


(15) 


1517.39631 


(0.00335) 


32(4) 


7.5 


0.34 




T 




658.6855 


(14) 


1518.17521 


(0.00322) 


33(4) 


7.7 


0.44 








658.2465 


(15) 


1519.18768 


(0.00341) 


31(4) 


7.2 


0.42 




5 




657.8284 


(06) 


1520.15327 


(0.00139) 


76 (4) 


17.8 






1 


* 


639.4081 


(06) 


1563.94636 


(0.00144) 


78(4) 


18.3 


0.48 






* 


638.9238 


(22) 


1565.13179 


(0.00543) 


21(4) 


4.9 


0.36 






* 


638.5688 


(11) 


1566.00196 


(0.00265) 


42 (4) 


9.8 


0.42 








638.1471 


(18) 


1567.03676 


(0.00446) 


25 (4) 


5.8 










603.9510 


(24) 


1655.76352 


(0.00670) 


19(4) 


4.4 




1 






595.5174 


(18) 


1679.21205 


(0.00494) 


26 (4) 


6.1 






6 




581.4725 


(25) 


1719.77176 


(0.00730) 


19(4) 


4.4 
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Table 4. Extracted pulsational frequencies, continued. 



Frequency 


Period 


Amplitude 


S/N 


Splitting 




n (= 2 






s 


ppm 










JJ7.7Z0D 


ZO) 


1785 94207 (0 00833 s ! 


17 C4") 


4.0 






t 


SS7 ASA/1 
JJ / .UJO'4 




1 793 ? 1 886 CO 00676 s ! 


22 (4) 


5.1 


1.38 




7 


SSA 77QS 




1 797 65640 CO 004 1 3 s ! 


36 (4) 


8.4 






1 






1 863 78608 TO 003 1 8 s ! 

IOUJ. / UVJV70 \\J .\J\J J i. 17 y 


50(4) 


11.7 


0.48 






S^A flASA 




1 865 44339 CO 00686 s ! 

10UJ.TTJJ7 ^VJ.UVJUOVJ,/ 


23 (4) 


5.4 








s 1 1 nrnA 

J 1 I.UUjt 


\\JI) 


r g^f. Q3497 co 00272 s ! 


64(4) 


15.0 






8 


/ISA /1Q7S 




2190 59091 CO 01022 s ! 


21 (4) 


4.9 








zl^S ^AQA 




2296 89908 CO 00970 s ! 


25 (4) 


5.8 








T.QO QA7S 

jyu.oo / j 


UZ) 


2558 41205 CO 00129 s ! 

Z,J JO.TlZ.UJ ^U.UVJIZ,'?^ 


233 (4) 


54.8 


0.29 




T 
1 


^Qfl S777 

J7U.J / Z / 


(ZU) 


2560 34269 CO 1 3 1 6 s ! 

V7 V/ ■ J 1 Z- V7 y VV7*V7±^/ XV7 / 


23 (4) 


5.4 






12 


^87 7^S 


UZ) 


2612 78405 CO 00124 s ! 

Z.VJ 1Z.. / UTuJ V"* 1 Z7-T y 


253 (4) 


59.5 


0.14 


T 

1 




"387 SQ1 ^ 
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J. H. Telting et al.: KIC 1 1558725: a 10 day beaming sdBV+WD binary 
Table 4. Extracted pulsational frequencies, continued. 
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Notes. The columns marked n t list the sequence number obtained from period spacing, for the t=\ and 1=2 sequences. Multiplet candidates used 
to derive the period spacings are indicated by the span of the arrows above and below the sequence numbers; for doublets only an up-arrow is 
used. The last column lists OA for orbital alias, or an asterisk if the frequency is part of a combination, as listed in Table 5. 
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J. H. Telting et al.: KIC 1 1558725: a 10 day beaming sdBV+WD binary 
Table 5. Possible frequency combinations 
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20 



Notes. All combinations F 3 =F 2 + Fi are significant to within the errors of the observed frequencies. The rightmost columns list the amplitudes. 
The top entry is an orbital alias (see Sect. 5. 1) 
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